All-optical diode action in asymmetric nonlinear photonic multilayers 
with perfect transmission resonances 
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Light propagation in asymmetric Kerr-nonlinear multilayers with perfect transmission resonances is 
theoretically investigated. It is found that hybrid Fabry-Perot/photonic-crystal structures of the type 
(BA)'^(AB)'^(AABB)'" exhibit both pronounced unidirectionality (due to strong spatial asymmetry of the res- 
onant mode) and high transmission (due to the existence of a perfect transmission resonance). This results in 
nonlinear optical diode action with low reflection losses without need for a pumping beam or input pulse mod- 
ulation. By slightly perturbing the perfect transmission resonance condition, the operating regime of the optical 
diode can be tuned, with a trade-off between minimizing the reflection losses and maximizing the frequency 
bandwidth where unidirectional transmission exists. Optical diode action is demonstrated in direct numerical 
simulation, showing >92% transmittance in one direction and about 22% in the other. The effect of perfect trans- 
mission resonance restoration induced by nonlinearity was observed analytically and numerically. The proposed 
geometry is shown to have advantages over previously reported designs based on photonic quasicrystals. 

PACS numbers: 42.65.Pc, 78.67.Pt, 42.25.Bs. 
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I. INTRODUCTION 

An optical device capable of unidirectional light transmis- 
sion (an optical diode) is a direct analogy to an electronic 
diode and has a great application potential in optical commu- 
nication and information processing. To this end, a variety of 
effects can be used, including magnetooptical WW , liquid- 
crystal p], nonlinear or gyroanisotropic ifioll elements 
in photonic microstructures. Also, several new concepts have 
been proposed recently such as indirect photonic interband 
transitions lEHl, planar chiral OIIl or negative-refraction 
|[l4ll metamaterials. Among all these approaches, nonlinear 
one-dimensional (ID) photonic multilayers deserve special at- 
tention because of their extreme simplicity in fabrication and 
modeling. This allows to explore the optical properties of non- 
periodic structures systematically and in a wide range of ge- 
ometries (see, e.g., 1 15]). In addition, a potential for versatile 
ultrafast pulse control by multilayer structures only a few mi- 
crons thick has been shown in recent experiments l fl6[[T7|] . 

The operation of a nonlinear photonic multilayer-based op- 
tical diode relies on two principles. First, the structure should 
be strongly nonlinear. A spectral resonance is usually em- 
ployed to enhance the nonlinear interaction, also providing 
the structure with a highly transmissive state which is very 
sensitive to parameter variations. Second, the structure must 
be spatially asymmetric so that the resonant mode couples 
more strongly to the wave impinging from one side (e.g., from 
the left) than from the other side. This causes a direction- 
dependent refractive index change due to the Kerr nonlinear- 
ity. So, the leftward-impinging wave can induce a sufficient 
shift in the resonant frequency cOres so that it matches the 
incident wave frequency cOq, resulting in a high left- to-right 
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transmission, while the rightward-impinging wave causes a 
weaker shift of (Or^s so that the (Ores 7^ resulting in a lower 
right-to-left transmission. Several design ideas have been pro- 
posed, achieving spatial asymmetry using graded refractive 
index H, 0] or asymmetrically placed defects (tI HI, [3] in a 
ID photonic crystal (PhC), or, more recently, using inherently 
asymmetric eigenstates in photonic quasicrystals [9] . Gener- 
alizations of the approach beyond the ID multilayer geometry 
have also been reported, based on quasi- ID coupled-resonator 
iilllil] or 2D ifllll PhC waveguides. 

Nevertheless, for a vast majority of photonic structures, 
the design principles for an optical diode seem to conflict on 
one important point. To increase the structure's sensitivity to 
the direction of incidence, one needs to increase the spatial 
asymmetry of the structure. However, to increase the maxi- 
mum transmission at a resonance peak, the structure should 
remain symmetric. It was once believed that only symmet- 
ric multilayers can exhibit theoretically perfect transmission 
resonances (PTRs) ll22l- [24]l an d that PTRs are explicitly re- 
lated to mirror symmetry [25] . Further studies revealed that 
perfect transmission is rare but possible in asymmetric multi- 
layers. Examples based on periodic i26ll . Fibonacci i27h . and 
Thue-Morse [28] geometry were given. In a series of recent 
works, Grigoriev and Biancalana [28, 29] report optical diode 
action based on PTRs in Thue-Morse multilayers. However, 
such structures typically need to consist of many layers due to 
relatively weak spatial asymmetry and field enhancement, in- 
creasing their sensitivity to the detrimental effects of material 
absorption and fabrication imperfections. Besides, quasiperi- 
odic structures possess a very rich variety of resonant modes, 
not all of which lend themselves to straightforward studies. 
Although there has been definite progress in analytical treat- 
ment of nonlinear resonances in Thue-Morse multilayers ll2^ . 
the overall picture remains complicated. 

In our recent work llsol] . it was confirmed that mirror 
symmetry is sufficient but by no means necessary to de- 
sign multilayers with PTRs. It was shown that structures 
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combining perfect transmission and highly asymmetric Hght 
locaHzation at the resonant modes can readily be achieved 
in a simple Fabry-Perot/PhC (FPPC) geometry of the type 
(BA)^(AB)^(AABB)^. Moreover, localization strength and 
asymmetry can be straightforwardly and independently con- 
trolled by changing k and m, respectively (see Fig. [TJ. It 
was envisioned that such structures would bring about an im- 
proved design of a nonlinear optical diode. However, an ex- 
plicit investigation of this design in the nonlinear regime has 
not been performed. 

In this paper, we investigate the nonlinear optical response 
of FPPC multilayers from the point of view of optical diode 
design. It is shown that FPPC structures are suitable for high- 
transmission and high-contrast unidirectional operation, re- 
sulting from combined contribution of PTRs and spatial asym- 
metry. A tradeoff between the strength of asymmetric re- 
sponse and the amount of reflection losses is established, so 
the best unidirectional operation is reported in a design where 
the PTR condition is slightly perturbed. Moreover, the effect 
of PTR restoration was observed where nonlinearity is seen 
to increase the maximum transmittance in a frustrated PTR 
back to unity. The optimized design exhibits optical diode ac- 
tion in the same input intensity range (around 10 MW/cm^) 
but with less than half as many layers as a Thue-Morse struc- 
ture reported in Refs. Hslllsl], owing to an increased localiza- 
tion strength in FPPC multilayers. By means of direct time- 
domain numerical simulations, optical diode action in the pas- 
sive regime was demonstrated with more than 92% transmis- 
sion. In the same parameter range, only pump-assisted or ac- 
tive operation was reported previously 1 18, 28, 29]. 

The paper is organized as follows. In Section |II1 the pre- 
requisites needed for PTR formation in FPPC structures are 
briefly reviewed. In Section |llll the nonlinear response of 
these structures is analyzed by using nonlinear transfer matrix 
method. The role of asymmetry and PTRs in optical diode ac- 
tion is identified. It is found that perturbation of the PTR con- 
dition leads to a tradeoff between reflection losses and unidi- 
rectionality, and that a proper structure design will cause PTR 
restoration in the nonlinear regime. In Section |IVl direct nu- 
merical demonstration of optical diode action in FPPCs using 
time-domain simulations is provided. Finally, Section|V] sum- 
marizes the paper. 



II. PERFECT TRANSMISSION RESONANCES 
IN MULTILAYERS 

We consider binary multilayer structures built of two kinds 
of dielectric layers (labeled A and B) with thicknesses Ja, 
Jb, linear refractive indices ^a, ^b, and Kerr nonlinear co- 
efficients Xa^ X^-> respectively, located in a homogeneous di- 
electric medium with n = hq. The building blocks A and B 
are assumed to be of the same optical thickness in the linear 
regime, so that 



riAdA riBdB Tic/ (2(0o) = Ao/4. 



(1) 



This equality is conventionally called the quarter- wave (QW) 
condition. It assures that the transmission T{(o) and reflection 
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Figure 1: (Color online) (a) Design and (b) transmission spectrum 

of an FPPC structure 8*^^'^^ as in Eq. (|4]), showing a PTR at cOq 

in the inset, (c) Asymmetric field distribution in structures S*^^'^^, 
S(3,5) 

(enhanced spatial asymmetry), and S*^^'^^ (enhanced localiza- 
tion strength) | 30]. The materials are chosen as in Ref. |29] to be 
polydiacetylene (9-BCMU, = 1.55, Xa = 2.5 x 10"^ cm^/MW) 
and rutile (Ti02, hb = 2.3, Xb ~ 0). The thicknesses and are 
chosen according to Eq. ([T]). 



R{co) spectra for an arbitrary arrangement of A and B layers is 
periodic in frequency with period 2(00. An additional mirror 
symmetry is present in each period, so that l l5L i3Q1 



T{co) = r((o + 2(0o), T{coo-Aco) = r((Oo + A(o), 
R{co)=R{co^2coo), R{coo-Aco)=R{coo^Aco). 



(2) 



Thus all the spectral properties of a given QW multilayer 
structure, determined by its layer arrangement, are contained 
in the region [0; cOq]. Note that QW multilayers always have a 
PTR at even multiples of (Oq, so that \T{2mcOo) \ — 1 exactly. 

When two arbitrary multilayers and 5*2 are stacked to- 
gether to form a combined structure S\S2, it is easy to relate 
its spectral properties to those of its constituents via Airy-like 
formulas at every frequency 



Ts,Rs2Ts^ 



Ts,Ts2 



^-Rs.Rs2 



(3) 



Here, Si denotes the inversion of 5"!, e.g.. Si = BAABA for 
5"! = ABAAB. Eqs. ([3]) can be used in a recurrent fashion 
with Fresnel formulas at dielectric interfaces acting as the ini- 
tial conditions, giving rise to an analytical way of calculating 
R{co) and T{co) for any multilayer structure [30]. 

It can be seen from Eqs. ^ that Rs^ {co) = Rs2{o)) = re- 
sults in Rs^s2{^) = ^nd 175-^5-2(^)1 = 1- That is, frequency- 
matched PTRs in the constituent structures always result in a 
PTR at the same frequency in the combined structure. It is 
known that both = (BA)^(AB)^ (a Fabry-Perot interferom- 
eter with a half-wave defect) and ^2 = (AABB)'^ (a ID PhC 
with doubled period) have a PTR exactly at CO = cOq. As a 
result, the combined FPPC structure 



■: S1S2 = (BA)*(AB)*(AABB)'' 



(4) 



also has a PTR at CD = fflo |30]. As the hght is strongly local- 
ized in Si (the Fabry-Perot part of the structure), the spatial 
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Figure 2: (Color online) (a) Linear transmission spectrum, (b) non- 
linear transmission spectrum around cOq showing PTR bending, and 
(c) input/output transmission hysteresis at ft) = 0.9935 ft)o for the 
structure 





0.9TO 0M5 1.000 i.O05 0.^90 0.995 1.000 1.005 0.990 0.995 1.000 

Frequency w/tjo 



86 0.993 i.OOO 

Frequency tL>/<U(j 





0.990 0.995 1.000 1.005 0.990 0.995 1.000 1.005 0.990 0.995 1.000 1.005 0.990 0.995 1.000 1.005 

Frequency io/ioq Frequency w/tJij Frequency cj/i(>o 



(b) Frequency 



Figure 3: (Color online) Nonlinear transmission spectra of (a) S*^^'^^ 
for = 5, 6, 7, 8 and (b) S^^''^) for m = 1 , 4, 7, 10. Red solid and blue 
dashed lines correspond to right-to-left and left-to-right propagation 
of light with input intensity 50 MW/cm^. Light gray upright peak 
is the linear transmission spectrum. Regions of prospective unidirec- 
tional transmission are highlighted as light yellow shaded areas. 



field distribution in the resonant mode is highly asymmetric. 
Moreover, larger k enhances the strength of field localization, 
w^hereas larger m increases its spatial asymmetry (Fig. [T]:), 
so they can both be controlled independently and straightfor- 
wardly. This makes the FPPC structures of Eq. ^ particularly 
attractive for a systematic study of unidirectional transmission 
involving a single-cavity resonant mode. In what follows, we 
are going to analyze the nonlinear optical properties of S*^"'^^ 
at its PTR around (Oq. 



III. NONLINEAR TRANSMISSION SPECTRA 
OF ASYMMETRIC MULTILAYERS 

To investigate the optical properties of FPPC multilayers in 
presence of Kerr nonlinearity, we assume that the refractive 
index varies with field intensity / as 

^fffcO = nA,B {l^XA,Bl{z^t)/2) (5) 

where / = I{z,t) varies both in space and in time, making 
the transmission spectrum intensity-dependent and potentially 
multistable 0, [isl HH] . To calculate it, we use the standard 
nonlinear transfer matrix method with layer subdivision | 311 
(note that more advanced methods have been introduced re- 
cently 1 29]). The example transmission spectra of S*^^'^^ are 
shown in Fig. |2l The PTR is seen to undergo characteristic 
resonance bending, in accordance with numerous earlier find- 
ings ill H [il 111 111, [11. Obviously, the resonance bends 
more strongly as the input intensity increases, and transmis- 
sion spectra eventually turn bistable (Fig. |2j)) with hysteresis 
behavior clearly visible on an output vs. input intensity dia- 
gram (Fig. 12:). Note that a slight dependence on the direction 
of incidence can already be noticed in the hysteresis loop, al- 
though in the effect is too weak and requires very high 
input intensities. 

From the FPPC design principles it is apparent how the 
changes in k and m should reflect themselves in nonlinear 



transmission spectra. Increasing k will increase both the res- 
onant mode 2-factor (making the resonance peak narrower) 
and the field enhancement factor (resulting in a stronger non- 
linear interaction and hence a stronger resonance bending for 
the same input field intensity). Increasing m should enhance 
the asymmetry in the field distribution, so that the nonlinear 
transmission of the structure becomes more dependent on the 
direction on incidence. This is confirmed in Fig.[3j For mod- 
erate values of input intensity (< 50MW/cm^) asymmetric 
transmission starts to manifest itself for ^ > 7, becoming more 
pronounced with increasing m. It can be seen that owing to 
direction-dependent asymmetry in resonance bending, there 
is a frequency range where a high-transmission state only ex- 
ists for one direction of incidence (see Fig.[3]3). In this range, 
high-quality unidirectional transmission is expected. 

Note in Fig. [3] that the maximum transmittance (Tiiax) in 
a bent PTR remains very close to unity despite the fact that 
Eq. ^ perturbs the QW condition O for ;^ 7^ 0. Hence 
the PTR condition (or, in fact, any recipe to provide a res- 
onance with near-unity transmittance) is important to ensure 
that the reflection losses of an optical diode will remain low. 
To demonstrate that the diode performance is very sensitive 
to the structure design, we compare 8*^^'^^ in Eq. (|4]) with its 
modification created by shifting the defect in its Fabry-Perot 
part, which can be expressed as 

(BA)^-^(AB)^+^(AABB)^. (6) 




1 = 2 



Frequency w/cjq 

Figure 4: (Color online) Nonlinear transmission spectra of the mod- 
ified structure as of Eq. ^ for k = l, m = l, and / = — 2 . . . + 2. Red 
solid and blue dashed lines correspond to to right-to-left and left-to- 
right propagation of light with input intensity 50 MW/cm^. 
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Figure 5: (Color online) Nonlinear transmission spectra of the mod- 
ified structure S'*^^'^^ as of Eq. for (a) positive and (b) negative 7 
ranging between -0.08 and 0.08. Red solid and blue dashed lines cor- 
respond to to right-to-left and left-to-right propagation of light with 
input intensity 50 MW/cm^. 



The results are shown in Fig. IH It is clearly seen that for 
the shifted defect (/ 7^ 0) the resonance is no longer a PTR 
and the linear transmittance comes nowhere near unity, re- 
sulting in poor Tmax in the nonlinear regime as well. On the 
other hand, it can be seen that the asymmetry in transmission 
gets drastically enhanced. So, while the optical diode action 
severely suffers in terms of reflection losses, unidirectional 
transmission becomes far easier to achieve. The reason is that 
increasing the number of periods in a Bragg mirror greatly 
increases its quality, in turn increasing the amount of time it 
takes the light to tunnel through such mirror into the Fabry- 
Perot defect mode. Hence, the difference of two periods in 
the (BA)^~^ (AB)^+^ structure not only destroys the PTR con- 
dition but also creates a much stronger asymmetric response 
than that introduced by the (AABB)^ part. This is under- 
standable because (AABB)^ is free of internal reflections at 
COo, in contrast to (BA)^^ 

Therefore, it can be concluded that FPPCs feature a trade- 
off between the maximum transmittance and unidirectionality. 
Hence it can be expected that an optimum performance would 
be achieved with a slight perturbation of the PTR condition. 
Consider a structure 



S/(/c,m) ^ (BA)^(AB)^(A'A'BB)" 



(7) 



where A' denote the A-layers with slightly changed optical 
thickness in the photonic crystal part, e.g., by 



^A' = (i + r^A. 



(8) 



For small 7, Fig. \5\ shows that the enhancement of unidirec- 
tionality more than makes up for a slight decrease in reso- 
nance quality. What is more, for 7 0.04 one can see that 
nonlinear Tmax becomes closer to unity compared to the linear 
regime, in the same way as it is seen for / = 1 in Fig.jH Hence, 
the S'^^'^^ design not only combines the benefits of structures 
given by Eqs. Q and © but also restores the PTR condition 
for a properly chosen 7. The reason is that the slight perturba- 
tion of the QW condition ([T]) by 7 is compensated for by the 
nonlinear refractive index shift ([5]) for a certain value of input 



field intensity. Fig. [6l comparing the 3D nonlinear transmis- 
sion spectra of 8'*^^'^^ vs. S^^'^\ illustrates that the asymmetry 
is much more pronounced for the former case, while Tmax re- 
mains close to unity in both structures. With the effect of PTR 
restoration taken into account, the S'*^^'^^ design is unambigu- 
ously superior in terms of the optical diode performance. 



IV. TIME-DOMAIN SIMULATIONS 
OF OPTICAL DIODE ACTION 

In this section, we demonstrate the unidirectional transmis- 
sion and optical diode action in optimized FPPC structures 
as given by Eq. (|7]) by means of direct time-domain numeri- 
cal simulations. We use the ID Lax-Wendroff method |[32ll as 
implemented in Refs. | 33, 34], which is analogous to a com- 
monly known Yee's finite-difference time-domain scheme. 

In accordance with the guidelines of the previous sections, 
S'(7,7) ^i^]^ 7= 0.05 is examined as a model structure. As seen 
in Fig.lSJ), this structure features a maximum linear transmis- 
sion of 91.2%. In the nonlinear regime, its right- to-left trans- 
mittance Tr increases towards unity due to the PTR restoration 
effect described above, while left-to-right transmittance 7] at 
the same frequency is low. The operating input intensity was 




Figure 6: (Color online) Nonlinear transmission spectra of (a) S*^^'^^ 
and (b) S^'^'^''^^ with 7= 0.05, for input intensities up to 55 MW/cm^. 
Two semi-transparent surfaces correspond to two different directions 
of light propagation. 
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with a rectangular envelope, long enough for the steady-state 
response to develop. The carrier frequency, determined by 
the input intensity values, was chosen at (0 = 0.9982(00 or 
12.1388 X 10^^ Hz. To avoid numerical artifacts caused by 
abrupt onsets of the rectangular envelope, Gaussian-like tran- 
sient excitations with decay time of 0.03 ps were applied im- 
mediately before turning on and after turning off the rectangu- 
lar quasi-monochromatic pulse (at 1 and 37 ps, respectively). 

The results are presented in Fig.[7t-b. It is seen that right- 
to-left transmittance 7> eventually reaches 92.05%, slightly 
exceeding the maximum linear transmittance of the structure 
(91.2%), although the transient processes are seen to last for 
up to 15 ps. In the reverse direction, the transmittance 7] re- 
mains below 23%. The on-off contrast can be improved by 
increasing both the operating intensity and (0, which will also 
decrease the transient time in the right-to-left propagation, but 
this will decrease Ty. The maximum nonlinear correction to 
the refractive index, monitored in the middle of the Fabry- 
Perot cavity during entire simulation, does not exceed the val- 
ues of 0.0085 (right-to-left propagation) and 0.003 (left-to- 
right propagation). These values are even lower than the esti- 
mation of 0.01 made in Ref. so the structure is expected 
to remain stable for the chosen intensities. 



Figure 7: (Color online) Results of FDTD simulations for pas- 
sive diode action in a S'*^^'^^ structure for /q = 9.5MW/cm^ and 
ft) = 0.9982ft)o: (a) with incident wave from the right and (b) with 
incident wave from the left. Solid lines, transmittance; dashed 
lines, reflectance. Blue (darker) color, left-to-right propagation; red 
(lighter) color, right-to-left propagation. In the background, quasi- 
monochromatic incident pulse is shown in light gray. The analyti- 
cally derived operating point of the diode in (c) intensity hysteresis 
loop and (d) nonlinear transmission spectra is also presented, in the 
same form as Figs.tSHU 

chosen at a rather moderate value of 9.5MW/cm^ in order to 
be able to work in the simple, passive regime. Higher input 
intensities (> 20MW/cm^), where the resonance peak gets 
sufficiently bent, make it increasingly difficult to excite the 
upper branch of the hysteresis loop by means of the incident 
wave alone. This can be overcome by using pump-assisted op- 
eration or modulated incident wave MlaSi], which IS 
rather unfavorable for optical diode applications but may be 
promising for all-optical modulator designs. 

The structure was discretized in the z direction with spa- 
tial mesh size Az = 5 nm, which corresponds to 30-50 cells 
per layer. The time step was determined by the condition 
that the Courant number V = 1, so = 0.01 667 fs. Such 
parameters ensure that the effects of numerical dispersion 
can be neglected in our ID calculations. Indeed, the test 
runs without nonlinearity confirm that the resonant frequency 
(1)0 = 27r X 12.161 X 10^^ Hz matches the results of linear 
transfer matrix calculations. 

Unidirectional transmission is demonstrated by examining 
the steady- state response of the nonlinear multilayer to a con- 
tinuous wave incident from either side of the structure. To 
simulate this, we use a quasi-monochromatic wave packet 



For comparison, the operating regime of the proposed opti- 
cal diode derived from the nonlinear transfer matrix method is 
presented in Fig.[7]:-d. Analytical calculation is seen to yield 
slightly better performance figures (7)- = 93.5%, 7] = 19.5%). 
These small discrepancies are attributed to spatial discretiza- 
tion present in both time-domain and nonlinear transfer ma- 
trix methods and can be seen as the effect of 1—2% struc- 
tural disorder. Other than that, analytical and numerical re- 
sults show a very good agreement. 

It is remarkable that the passive unidirectional transmission 
is observed at the lower-intensity "crest" of the hysteresis loop 
(see Fig.[7J:) rather than at its higher-intensity "tail" as in ear- 
lier works [18]. Previously, only active or purnp-assisted op- 
tical diode action was achieved in this regime |[2Tl[28[ [29l]. 

Note, finally, that the demonstrated optical diode action oc- 
curs in the same input intensity range as in the recently re- 
ported Thue-Morse multilayers [S, '29] but with the use of 
a structure less than half as thin and containing half as many 
layers (56 versus 128). This results from an increased field 
localization strength in FPPC structures. Another advantage 
of the FPPC design is that it allows to purposely engineer a 
PTR-featuring resonant mode with desired properties, rather 
than rely on naturally occurring PTRs in photonic quasicrys- 
tals where not all resonant modes are suitable for unidirec- 
tional operation (see |(29!]). The operating intensities can be 
brought further down by increasing k, with a relatively in- 
significant incease in the number of layers. However, this 
will also increase the length of the transient response of the 
structure, effectively resulting in a tradeoff between energy 
efficiency and maximum attainable bandwidth when such an 
optical diode is used in communication systems. 
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V. CONCLUSIONS AND OUTLOOK 

In summary, we have theoretically investigated the optical 
properties of asymmetric Kerr-nonlinear multilayers with per- 
fect transmission resonances. Using the nonlinear transfer 
matrix method, it is found that FPPC structures of the type 
(BA)^(AB)^(AABB)^ exhibit both pronounced unidirection- 
ality and high transmission, which makes these structures suit- 
able for nonlinear optical diode action in the passive regime 
with low reflection losses. A trade-off is found between the 
maximum transmission and unidirectionality, subject to tun- 
ing by slightly perturbing the PTR condition. An effect of 
PTR restoration is observed in perturbed FPPC structures. 
Theoretical predictions are confirmed in direct time-domain 
numerical simulations, yielding more than 92% transmission 
for the input intensity of <10 MW/cm^. This value exceeds 
the linear maximum transmittance of 91.2%, further confirm- 
ing the PTR restoration effect. Passive optical diode regime 
was demonstrated where only pump-assisted or active opera- 
tion was reported previously JTH |28l |29|], and the proposed 
geometry offers the same performance as recently reported 
Thue-Morse multilayers [28^ 29J but with less than half as 



many layers. 

The results obtained can also be viewed as a systematic 
investigation of a single photonic-multilayer resonant cav- 
ity in the nonlinear regime. They can be used as a starting 
point for investigating more sophisticated multilayer geome- 
tries combining spatial asymmetry and high transmission. For 
example, recent studies report that coupled-cavity resonances 
in multiple-defect PhCs |35] and photonic quasicrystals |29J 
have very different character of nonlinear bending and dif- 
ferent properties of unidirectional transmission, compared to 
single-cavity resonances. To examine PTRs created by mul- 
tiple Fabry-Perot elements coupled in an asymmetric fashion 
would therefore be a very interesting extension of the princi- 
ples derived in the present work. 
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